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Abstract 

The generalization of the Yang-Baxter equations (YBE) in the presence of Z 2 grading 
along both chain and time directions is presented and an integrable model of t — J type 
with staggered disposition along a chain of shifts of the spectral parameter is constructed. 
The Hamiltonian of the model is computed in fermionic formulation. It involves three 
neighbour site interactions and therefore can be considered as a zig-zag ladder model. The 
Algebraic Bethe Ansatz technique is applied and the eigenstates, along with eigenvalues 
of the transfer matrix of the model are found. In the thermodynamic limit, the lowest 
energy of the model is formed by the quarter filling of the states by fermions instead of 
usual half filling. 



LAPTH-804/00 
hep-th/0006243 
June 2000 

1 e- mail : am bjorn @alf. n bi . dk 

2 e- mail : arnaudon @lapp . in2p3 . fr 

3 e-m&i\:sedrak@lx2.yerphi.am, Permanent address: Yerevan Physics Institute, Armenia 
4 e-m&i\:tigrans@moon.yerphi.am, Permanent address: Yerevan Physics Institute, Armenia 
5 e-mail : sorba@lapp.in2p3. fr 



1 Introduction 



The interest to ladder type models was raised in a beginning of 90-s (see for a review 
lUl) in connection with high temperature superconductivity problems in metal oxides. It is 
believed that quasi-one dimensional multi-ladder chains of strongly interacting electrons 
reflects the most important aspects of two dimensional systems and also can reveal some 
properties of the week coupling between conducting planes. 

Recently there has been considerable interest in the construction of integrable ladder 
type models motivated by the desire to use powerful technique of Algebraic Bethe Ansatz 
(ABA) 0, |3|] in the exact investigations of the variety of physical phases of the models. 

In the articles || integrable ladder models were constructed by extension of symmetry 
algebra, in |J by defining first the ground state and then formulating a model, which has 
it. The higher conservation laws of an integrable models, which contains next to nearest 
neighbour interactions, were used in construction of ladder models in JF| by developing the 
approach of the article ||. There are also some other type attempts in this area || |I0| . 

Usually integrable models are homogeneous along a chain, namely, the spectral u 
and model parameters are the same in the product of i?-matrices along a chain. It is 
obvious, that if one considers arbitrary shifts of the spectral parameters by some zi in the 
monodromy matrix we still will have an integrable model. But in order to have a local 
Hamiltonian we need to consider shifts with the fixed periodicity n, which will cause the 
interaction of spins (or electrons) within an amount of n neighbours leading to n-ladder 
model. This type of model was first considered in [11], developing the ideas of the articles 
|T2|1 (see also then in a chain of articles [JT4J, [L^, [Uj. 



In the article [|J we have proposed an inhomogeneous model based on XX Z spin-chain, 
where inhomogeneouty appeared not only in the staggered shifts of the spectral parameter, 
but also in change of structure of .R-matrices in a product along a chain. Namely, two 
monodromy matrices of chains M s , s = 0, 1 were considered along a time direction, where 
the i?-matrices in the product have an alternating disposition of the anisotropy parameter 
±A of the XX Z model. Besides that, and contrary to case considered in |T^, TJ, H), HJ, 



the spectral parameter of the second line has an opposite sign. Due to the double space 
translational invariance the Hamiltonian of the model contains interaction between three 
neighbours sites of the chain and therefore represents a zig-zag type ladder model. At the 
free fermionic point A = the model becomes a model of two noninteracting fermions, 
hopping separately in the odd and even sites of the chain. 

In this article we are extending our approach to spin — 1 case and considering the t- J 
model [O, BTL p3l |24|, ER |26[ . Following H we are writing two Yang-Baxter equations 



(YBE) (see equations || and below) for each step of two alternating /^-matrices in the 
chain, but as it appeared in the solution, in this isotropic in the spin space model two 
intertwining i?-matrices are coinciding. As a result we are being left with the alternating 
shift of the spectral parameter and the change of sign of the spectral parameter in second 
line. As we will see, though this gives us the same Bethe Equations (BE) for the spectral 
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parameters of the excitations as for the model derived in [[TJJ, 13, [TE| , their energy and 
the energy of ground state are different. It appears that in the thermodynamic limit the 
lowest energy of the model is reached by the quarter filling of the states by fermions, 
instead of the usual half filling. 

In Section 3 we formulate the model and find a local Hamiltonian, which has a zig-zag 
ladder form. It consists of two chains with t — J type Hamiltonians on each of them, the 
hopping term of electrons from one chain to the other and two type of interaction terms 
between chains. The first interaction term has the form of spin-spin interaction, where 
one spin is composed from two fermions on the same site of the chain, while the other 
spin is composed from two different fermions on the neighbour sites of the other chain of 
the ladder. The second interaction term has a topological form of interacting spins and 
is written for the triangles consisting of the zig-zag rungs. 

In Section 4 we apply ABA in order to find the eigenvalues and the eigenstates of 
the model. At the end of Section we find the ground state energy of the model and the 
spectrum of excitations in thermodynamic limit. 



2 The Yang-Baxter Equations and its solution. 

Following [f|] let us define the monodromy operator of the model as a product of the 
two chain monodromy operators M (u) and Mi(u) 

M(u) = M Q {u)M x {u), (1) 

where 

Mx(u) = Rt{u)R^ li2 {u)R^{u).... 

M {u) = R 01 {u)R^{u)Ro 3 {u).... (2) 

and we have two operations L\ and l 2 , defined with property i\ 2 = i 2 2 = 1. 

The transfer matrix of the model will be defined as a trace over auxiliary spaces and 

0' 

T (u) = Tr Tr 'M(u) = t (m)ti(m). (3) 

As it was shown in the article |4| the commutativity of transfer matrices t(u), t(v) for 
arbitrary values of the spectral parameters u and v can be ensured by the following set 
of Yang-Baxter Equations 

R 12 {u, v)R 2 \{u)R l2 {v) = R L 2 %{v)R 12 {u)R 23 {u, v), (4) 

and 

R 12 (u,v)^(u)R^ 2 (v) = R^\v)R? 2 (u)R 23 (u,v), (5) 
where we have used the braid (check) formalism for the convenience. 
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It is also convenient to use the fermionic operator formalism for R-matrices and the 



YBE developed in [18, 19, 20, 17 



For the t — J model under consideration |T9| , pl| , [23| , |24| , |25| , |26| we have to realize 
the spin 1 algebra in the fermionic approach. A minimum of two sorts of fermions is 
needed in order to express three basic states | +), | 0), | — ) of the spin — 1 particle with 
the z component of the spin equal to 1, 0, —1 correspondingly. 

Now let's define c a + , c a , where o =tl, as a creation-annihilation operators of fermions 
with the up and down spins respectively, together with their Fock space | 0),| a). 

The states with definite third projection of the algebra SU(2) can be realized through 
fermionic Fock space as follows 



-)=|0,|),|+)=|t,0),|0)=|0,0), 



(6) 



numerated as | 1), | 2), | 3) respectively. 

As it is obvious from the formulas || we have constructed a graded space with the 
following parities for the basic vectors 



p(I+»=p(I-» = i, P(|0» = 0. 



(7) 



In order to proceed further and write the fermionic .R-matrix we should calculate the 
Hubbard operator X™ =\ m)(n |; m, n — 1, 2,3 





/ i-x-i 


l-X+l 


l-Xoi 




l+X-l 


l+X+l 


l+)(0| 




V lo)(-l 


|0)(+| 


|0)(0| 



(1 - n^)n [ c+c T (1 - n T )c+ 

c|q n T (l-n|) c+(l - n{) 

(l-n T )q c T (l-n|) (l-n T )(l-n|) 



The trace of this operator is 



A = X™ = l-n T n i , 



(9) 



which is an identity operator on the space of states,where the double occupancy of the 
sites by fermions is excluded. 

Following [|]1J let's write down the fermionic i?-operator for the t-J model {spin 1 



I] model). 

R id {u) = a(u)I hJ + b(u)U hJ = a(u)I hJ + b{u) £ {-If^X^X^ 



N 



(10) 



m,n=l 



where H^- is the graded permutation operator of the spaces Vi and Vj. 
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Now by putting the i?-matrix form flTU|) into the YBE (H) and ([5]) and after some 
calculations one can find 12 equations which require the operations ~ and transformation 
ii to be 



a(u, v) = a(u, v), b(u,v) = b(u,v), 
a i2 (u) a(u) 



'ID 



b L2 (u) b(u)' 

The conditions ([II]) are reducing the 12 equations to following two equations 

a(u, v)[a(u)b(v ) — b(u)a(v )] + b(u, v )a(u)a(v ) = 
a{u, v) [a' 1 (u)b Ll (v) - b Ll (u)a Ll (v)} + b(u, v)a Ll (u)a Ll (v) = (12) 

the consistency condition of which can be found easily as follows 
b(v) b Ll (v) b(u) b Ll (u) 



a(v) a H (v) a{u) a il (u) 



constant = 6. (13) 



Here 9 is the new parameter of our model. 

Then the solution of fll2|) for the intertwinner parameters a(u, v) and b(u, v) will be 

Ku,v) = 6(u) _ b{v)_ ^ 
a(u,v) a(u) a(v) 

We should now define the i\ operation. It is easy to see from the expression ([13]) for 
v = u that i\ operation can be consistently defined as follows 

b L1 (u) _ b(u) b Ll u _ b{u) ^ 
a Ll {u) a{u) ' a Ll (u) a(u) 

It is clear from the formulas (0-0) that the ratio can be taken as a spectral 
parameter ^4 = u, = u = 9 — u. 

1 a(u) ' a(u) 

Finally, after appropriate normalization of a(u) and b(u) in order to have a(u) + b(u) = 
1, one finds the following solution of YBE 

a(u) = , a(u) = , a(u,v) 



1+u y ' l+9-u y ' 1+6-u-v 

, u -. . 9 — u ,. , 9 — u — v , , 

K„) = _, Hu) = TTJ - u K*,v)= 1 + 9 _ u _ v (16) 

According to standard prescription of the ABA technique the logarithmic derivative 
of the transfer matrix at some point defines the Hamiltonian of the model 



As it is known, in order Hamiltonian to be local, it is necessary to have a value uq, 
such that 

/.',..; ('/„) /,,- (18) 

Analyzing the solutions of our YBE ([H]), one can see from (P~0| ) that, at the point uq = 0, 
only Rij(0) = and Rij(0) = Iij. As calculations shows the Hamiltonian is nevertheless 
local, but it contains interaction between four neighbour points. 

Technically, in order to calculate the logarithmic derivative (|17D , one should put the 
expression of Rij operators around uq = up to linear terms 

Rij(u) = Ri,j(9) - uH id , 
Rlj(u) = Iij - uHij, 
R*(u)'=Rij(-9)+uH itj , (19) 

with 

m,n 

into the expression (|TJ) of the Monodromy Matrix 

M(U) = R^ (u) R\2 {u) ^23^34 («)••• -R\2* («) #23 (u) R^ 2 (u) #45 . . . . (21) 

As a result, after some algebraic calculations, we will obtain the following Hamiltonian 
for the present staggered t-J model 

JV 



H = 0A E E( 2 -y-y) " <£w) 

i=l U=U 

+ E[( 1 - ! T i + (-!)^) <Wv - (l - ^ - (-1)^) ctc i+2 , a 
+ 2 [(S 

+ 0(-l) 4 fei-iSw - \n*-m l+1 + w *-i+^+^ _ ie ^StS h i+1 S c i+2 } A, 

where A = Y\!i=\ A?; = l~Iili(l ~~ n ^ n i[) is the projector which excludes the double occu- 
pancy by electrons at any site i. 

The spin operators Si and S it j are defined as follows 

Si = = \^4,ji4^ 

Sij = \*tjtfW = i^c&d^Aj, (23) 
5 



where 

^ = (l-n T )q, * 2 = c T (l-7i i ), (24) 

and a are Pauli matrices. 

3 Algebraic Bethe Ansatz for the staggered t-J model. 



In this section we will apply the technique of ABA [BL |25|, |26j to the present model 



and find the eigenvalues and eigenstates of the Hamiltonian (|22|) . 
For this purpose let's introduce the L operators as follows 

(Li,j)k = (k I Rij I k') (25) 

which is a matrix in the horizontal auxiliary space and an operator in the quantum space. 
In matrix form it looks like 

a(u) — b(u)(l — n-f)nj, — b{u)c^c^ b{u){l — n^)ct 

U,j = [ -6(tx)cJ"c; a{u) - 6(u)n T (l - n;) &(u)cf (1 - n 4 ) 

6(u)(l-n T )q. &(w)c T (l-n|) o(u) + 6(u)(l - n T )(l - nj 



(26) 



The monodromy matrix Mj°,(u), which defined by matrix elements of the monodromy 
operators ([L],|2]) in the auxiliary space can be expressed as a product of L it j matrices as 
follows 

M Q {u)i = (k i m (u) i k<) = (-i)^^(i&){ 1 (AB)£....(A W )ir- 1 , 

Mx(«)* = (A; | Mi(w) | k') = (-lr^'H^)^^)^---^)^ 1 - (27) 
Following the notations of the article PEL one can write 




A S) u(u) A S) i 2 {u) 
M s (u) k k , = | A s , 12 (u) A s , 22 (u) B s , 2 (u) | s 0. I (2t 
C sA (u) C Sj2 (u) 

where A StCL b, B s>a , C g)Cl , -D s ; (a, b = 1, 2) act on the quantum space. 

Then, as a super-trace of the Monodromy matrix (]28D the transfer matrix (§) will have 
the form 

r s (u) = -A, i ii(u)-A S)2 2(u) + D a (u), s = 0,l. (29) 
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This matrix elements of the Monodromy matrix are obeying the algebraic relations 

= (_l)rfV>(*mH*(mO)^ 

( _ 1)pOTWm ') + K m "))^ m/ (u,v)M™i„(u)M( k „(v) 

= (-l) p(fc ' ){p(m)+p{m,)) M >( W )Mf ife; (n)^,,(n, v), (30) 



which are the consequence of YBE (|}||). In getting fl3~0|) we have used the properties 

(O). 

Consider now the empty fermionic state as a test "vacuum" 

| n) a =\ 0,0,...,0) s =| 0) u I 0) 2s ... I 0) Ns , s = 0,l, (31) 

and let's check that | Q) indeed is a eigenstate of transfer matrix (p9|) 

t s (u) | fi) a = | ft>!_ s (32) 

From ( P7|) t s (u) is a product of L it j matrices. Hence, in order to check ( |3"2] ) we should 
first calculate L fc | 0) s ,k- It appears that 

/ b(u)^ & t2 («)4| \ 

^|0> S)t = b*(u) b*(u)c+ s = 0,l. (33) 

\ a t2 (w) + b L2 (u) ) 

We see that | 0)^ is a upper-triangular matrix. Therefore the action of the product 
of Lofc in the formula ( p7|) on vacuum | f2)& as a matrix will also have upper triangular 
form 

[b^( u )b' 2 {u)]i B 1A (u) 

1 

[b^(u)b{u)]f S ,i(w) 
Mo(u)^|fi) = ( [b^(u)b(u)]f B 0i2 (u) ]|fi) , (34) 

1 

where we have used that a{u) + 5(w) = 1. 

We see that the B s i{u) and B s2 (u),(s = 0,1) operators create one particle states 
while C Sj i(u), C S fl{u) operators annihilate them 

C s , a (u) \Q) S = 0, s = 0,1; a = l,2 (35) 

We see from the expression (^4[) that 

v V)(u) = l-2[b tl (u)b(u)]% , (36) 
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where b Ll (u) and b(u) denned by the equations ( ]TB| ) and N is the length of the chain. 
This observation leads us to the following Ansatz for the eigenstates of t(v) 

| v 1 ,v 2 ,...v n I F) = F a -- a 'B 0>ai (v 1 )B lm (v 2 )...B 0>an (v n ) \ Q) , a % = 1,2; (37) 

is a n particle state. The F an ~ ai is a function of spectral parameters Vj to be specified 
later. 

The action of the transfer matrix fl29|) on the states fl3T|) is determined by the relations 
(0) and the intertwining properties of the A Syab (u), D s (u), B s>a (u) operators defined from 
the YBE (0). The components of the intertwining relations, which we need for the 
construction of the ABA are 

Di(u)B 0a (v) = ^—B 1>a (v)D (u)-^^B 1>a (u)D (v), 
b{u,v) b{v,u) 

A 1M (u)B 0tC (v) = % ^ B 1>c ,(v)A 0>b , a (u) + ^i^B 1>b (u)A , ca (v), 

o{u, v) b{u, v) 

B 1>a {u)B 0tb (v) = ri a '{u,v)B lta/ B 0tb ,{u), (38) 

where 

rl c \v) = -a(v)5i5 b ; + b{v)5«5i = -a(v)I^ b ' - b(v)U { b f c ' . (39) 

is a graded p( 
fulfills the following YBE 



Here n£ ),6 ' c is a graded permutation operator for p(l) = p(2) = 1, one can check that it 



r(A-/023r(A)££r(/i)*5 = r^)%%r{X)^r{X-^X- (40) 

Now by use of (PH), we can obtain, that the diagonal elements of the monodromy matrix 
act on the states ( p7|) as follows 

1 

Di(tt) I v u v n I F)i = TT ~i V I u i> u ™ I F )o + 

n n 

+E( A *)t:>w«) n ^(«i)i«>o, 

fc=l j = ljj£k 

[A 1>u (u) + Ai i22 («)] I «i, w„ | F)i = 

n 1 n 
01 It, Vi) 

1=1 v ' 1=1 

n n 

+j2( A ^± Fan --- aiB ^M n B >M)\n)o, (4i) 

fc=l i=l,jy^k 
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where 

= str[Z n (ti,v n )Z(n- l)(w, v nr .i)...l 1 (u 1 «i)], (42) 

and 

P*(«, = r J^« fc («' V *)- ( 43 ) 

As it follows from is a 2 x 2 matrix, the elements of which are the operators 



/ ( ni \ _ f Z M l li \ - f b ( u ) 1 ~ a ( u ) e l -a{u)e\ \ . . 

where e b a are quantum operators in the n-th space with matrix representation (e^)a = 5° d>L 
It is obvious, that the eigenvalue condition 

{D s>1 (u) - A S;ll (u) - A St22 (u)) | v u ...v n | F) s = u s (u,v 1 , ...,v n ) | vj, ...,t>„ | F)i_ s (45) 

will be fulfilled if 

i) we impose the cancellation of unwanted terms in fl4"T|) 



[(A fc )ti - Wa\:: b z] Fan - ai = ( 46 ) 

called Bethe equations (BE), and 

ii) we solve the eigenvalue problem for the small transfer matrix fl42|) 



r ii )a k' a "( M ;^ -,v n )F ai - an = uW(u,Vi)F<-<, (47) 
then we have the following expression for eigenvalues 

n ^ " 1 

^(u;«i, ...,«„) = (I- r + [6' lt2 (n)^( M )]^ J]- rvVfaVi), 



n 

vq{u;vi, ...,v, ' 



TT— ^— + [b t2 (u)b t2 (u)]%f\ -, 1 V 1 ^,^). (4f 



For the solution of second equation ii) we should make the ABA for a small auxil- 
iary problem of chain, with length n (number of particles) and "nested" transfer matrix 

T^. a .a n '""(«; Vi, v n ). This is why all this procedure is called Nested Algebraic Bethe 
Ansatz (NABA). In the article P5| |, it was demonstrated how to calculate A**'"^ and 



-^■aV.'a" an d to reduce the condition of cancellation of the unwanted terms for the ordinary 
t — J model to some equation. It is not necessary to repeat the same calculation here 
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since it differs very little from the carried one. The only difference is appearing in the 
term 

(A ltll + A 1>22 ) | fl> = [&«»P 2 («)]£ | Q), (49) 
therefore we obtain the following conditions 



Bethe equations. 

In the next step of the NAB A we should find the eigenvalues and eigenstates of t^ 1 ' (u). 
It is clear from the equations ( pEUD that we have another small integrable model with the R 
matrix [u) defined by the formula (|39|) and the corresponding transfer matrix t^ x >{u). 

Therefore we should apply a non ordinary ABA to this problem. The Y BE for the 
problem is 

rtf(u,v)MP a ' \u)M^ \v) = M^ b \v)M^'(u)r^''(u,v) (51) 

where Ma^ a is the corresponding (nested) Monodromy matrix. 
Now if we define 

M(1) ^ = ( d^U ) > r(1) ^ = - A(1) ^ - D(1) ^ w 

then by use of the formula ( |3"9"D and YBE (pT|), we find 

dM(u)bU(v) = -^B^\v)D^\u) + ^^B^(u)D^(v), 
b(u, v) b(v, u) 

A^(u)B^(v) = ^^B^(u)A^(v) + -^—B^(v)A^(u), 
b[u,v) b[v,u) 

BV)(u)bV>(v) = BV>(v)BV>(u). (53) 

Let's take as reference state 

I o>« = ( J 

n 

l^) (1) =|0)l 1) ...|0)« = (9)|0)l 1) . (54) 



k=l 

The action of the nested monodromy matrix M^(u) on the reference state | Vl)^ is 
described by the action /&(«) on | 0)£\ which we can find from (|44|). So we obtain 

n 11 hi \ 

a^\u) i n>« = - acu,^)] i = n ir^H i ^> (i) > 



i=l 



d {i \u) i = n & K^o i ( 5s ) 
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For the eigenstates of t^'(v), we are going to do the following Ansatz 

| «« . . . , «£>) = flW(^)BW(^) . . . S«( V £)) I (56) 

The action of t^ 1 ) (it) on the states (|56|) is the same as the action of the diagonal elements 
of (B2p on that states. By use of (RSI) we will obtain 



m m 



k=l 
m 



am(u)\ v?,... ,„£)> = n 



1 fr&M, (i) 



n 



(57) 



fe=i 



From the expression (|57|) we can easily write the eigenvalues of r* 1 * 1 ( 



u 



fo(it, 



m ^ n 

n^r^n feK . )U) 

m ^ n 

+ n 7-^1)711^) 



=1 6(«> v i ) j 



i=l 



y l ) ■ ■ ■ ) y m / ■ 



(5? 



One can get simply the first set of Bethe equations by comparing (|58| ) with the formula 
j). Inputing u = v & in (|50f) we obtain 



[6 11 («fc)5(«*)]* = n6(t;i 1) ,t; Jfc ), A; = 1,2,. 



n. 



(59) 



i=l 



The second set of Bethe equations, which are the conditions of cancellation of the 
unwanted terms and are similar to the corresponding equations of the standard 
XXX model and can be found easily as 



3=1 



; * J-ll, r (i) ay k-L,2,...,m. 



(60) 



This is exactly the same equation as found in [26 
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Finally we find 



{Vi}) 



m ^ n 

II w — mrll 6 ^) 



n 



\ b(vi,u) 



- [b Ll {u)b(u)) 



n 



n 



m ^ n 



(61) 



as the n particle state eigenvalues of transfer matrices Ti(u) and tq(u) respectively. 

But the transfer matrix of our staggered model is the product of tq(u) and Ti(u) 
therefore the eigenvalues u(u, {vi}) of r(u) are 



u(u, {Vi}) = u (u, {Vi})vi{u, {Vi}) 



(62) 



with the Bethe equations (|59| ) and ( j60| ) unchanged. 

Let us now to calculate the energy of excitations over the test "vacuum" | fl), called 
bare energy, which will be dressed in a real ground state due to interactions with particles 
in a filled Dirac sea. The bare energy is a logarithmic derivative of eigenvalues (|6l|) and 
( p2f) at the point u = 0. The simple calculation gives the energy and the momentum of 
n-particle state | vi, ...v n \ F) as it follows 



£<■({«;}) = "E 



t»| + 1/4 ( Vj - ey + 1/4 

+ 1/2 
0-1/2 



(63) 



where we have redefined the spectral parameters as u,- —*Vj — 1/2. 

The solution of the (|59D and (|60D is usually obtained in the thermodynamic limit 
(N, n, m — > oo, with the fixed ratio ^). In this case instead of a discrete set of spectral 
parameters Vj one introduces the distribution of continuous density p{v) of rapidities. The 
ground state is defined by filling up the Dirac sea(s) of negative energies by the electrons. 
It was argued in the article [25] that the ground state of the t — J model is defined by 
the string solutions of length two, which are filling of all states with negative energy. The 
lowest energy value can be reached by maximal filling of negative energy states, which 
corresponds to = 1 and with zero magnetization, corresponding to m — |. In our 
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model it is clear from the expression of the energy ( |6"3"D that only the spectral parameters 
of the interval 

— oo < u < - (64) 

have to be filled in order to form a ground state. But since this is exactly equal to half of 
lattice sites N, we will have a ground state corresponding to quarter filling of all states. 
After introducing the densities, the BE becomes an integral equation of the form 

Mv) + I Q du (v- Pi u)Ki = I {^TITi + ^-^ + 1/4 } ' (65) 

where Q defines the rapidity of the Fermi level. In case of quarter filling it is equal to 
zero. 

The energy of ground state is defined by the equation 



f 00 

En = -2N 



1 



v 2 + 1 (v- O f + 1 
= -2N(p(0)-p(6)). (66) 

At usual half filling Q = — 00, and as it follows from the expression (|66D , E = 0. 

It seams to us, that this model provides an interesting possibility to analyze by means 
of exact integrability the physics of systems with the quarter filled ground state. 
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